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Abstract 

A sixth order quadrupole boson Hamiltonian is treated through a time dependent variational 
principle approach choosing as trial function a coherent state with respect to zeroth 6q and second 
b\ + h^_2 components of the quadrupole bosons. The coefficients involved in the model Hamilto- 
nian are chosen so that the classical effective potential energy term has two distinct minima. The 
classical system is described by two degrees of freedom and has two constants of motion. The 
equation of motion for the radial coordinate is analytically solved and the resulting trajectories are 
extensively studied. One distinguishes three energy regions exhibiting diff'erent types of trajecto- 
ries. When one passes from the region characterized by two wells to the region of energies higher 
than the maximum value of the effective potential the trajectories period exhibits a singularity 
which reflects in fact a phase transition. The classical trajectories are quantized by a constraint of 
the integral action similar to the well known Bohr-Sommerfeld quantization condition. The semi- 
classical spectra corresponding to the two potential wells have specific properties. The tunneling 
process through the potential barrier is also studied. It is a remarkable fact that the transmission 
coefficients exhibit jumps in magnitude when the angular momentum acquires certain values. One 
suggests that for such angular momenta the rates of alpha or heavy cluster emissions might be 
very large. 
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I. INTRODUCTION 



Description of spectroscopic properties for a many body nuclear system is one of the most 
interesting aims of theoretical nuclear studies. Some of such properties are obtainable by 
using only few degrees of freedom. Thus the liquid drop model |l| tries to describe the basic 
quadrupole collective properties by using a harmonic liquid drop whose surface vibrates 
with small amplitudes around an equilibrium spherical shape. Few years later another two 



solvable simple formalisms have been launched namely the Jean-Wilets model for gamma 

an 
and the Davydov Filippov (DV) model of triaxial rigid rotor 3]. Although 

these models are not able to describe the complex details of nuclear spectra they have the 

big merit to define the framework and to place reference landmarks in the field. Indeed, 

many of the improvements proposed interpret their results by considering the departure 

from the reference picture provided by the simple models mentioned before. Thus, the first 

phenomenological model which goes beyond both the liquid drop and DF models is the 

vibration rotation model jj] where the drop surface may oscillate around a static deformed 

shape. The anharmonic terms in the quadrupole shape coordinates are added to the liquid 

drop for the first time by Greiner and Gneuss ^. This anharmonic model was extensively 

used to describe various types of spectra , gamma unstable, gamma stable, triaxial shapes 

in several publications [6]. The drawback of this model is the large number of parameters 

involved although they depend smoothly on the atomic mass number. A much smaller 

number of parameters is used by the coherent state model (CSM)^] which treats three 

interacting rotational bands within a restricted collective space defined through angular 

momentum projection technique, from three deformed states: one of them is describing 

the deformed ground state and is of an axial symmetric coherent function and the other 

two are simple polynomial excitations of the ground state determined by requiring the three 

states to be mutually orthogonal. GSM exploits the semi-classical features of coherent states 

which results in obtaining a suitable description for states of high angular momentum. It 

also provides a unified description of near vibrational, transitional and of well deformed 

nuclei. Another interesting phenomenological model is the interacting boson approximation 

which uses besides the quadrupole bosons a monopole one The model Hamiltonian is 

conserving the total number of bosons which as a matter of fact is limited to half the number 

of valence nucleons from the nucleus under consideration. By construction, the Hamiltonian 
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accounts for one of the symmetries 0(6), SU(5) and SU(3). Consequently, the nuclear states 
are just the irreducible representations for these extreme symmetries. However, to describe 
nuclei departing from these symmetries additional terms are to be added which blurs the 
primary beauty of the model. 

The phenomenological descriptions involving anharmonic boson Hamiltonians make use 
either of approximate methods or of diagonalization procedures to get the desired eigen- 
values. In the latter case the procedure is difficult to be applied when we are close to the 
picture where the results are very slowly convergent against increasing the basis dimension. 
On the other hand, when an approximate approach is used one does not know how reliable 
it is . In other words, we don't know how far the obtained result is from the exact one. 

The same difficulties are encountered also in microscopic formalisms. The microscopic 
theories interpret the experimental data in terms of single particle motion. Thus, the collec- 
tive motion is defined in terms of individual degrees of freedom. The random phase approx- 
imation (RPA) or boson expansions are defining boson operators in terms of bifermionic 
operators and then the Hamiltonian and transition operators are written as series of the 
RPA bosons ^, ^ . The convergence properties of the boson series and the magnitude 
of the ignored "rest", when a truncation is performed, can be evaluated only for ex actly 



solvable models. The most known solvable microscopic models are those of Moszkowski 12| , 
Lipkin-Meshkov and one level pairing Hamiltonian Q]. Applying for example, one of 
the above mentioned approaches to a solvable model and comparing the results with the 
corresponding exact output one may conclude upon the vahdity of the procedure adopted. 
It should be stressed the fact that using a realistic Hamiltonian instead of a solvable one, 
one expects that many of the features are at least qualitatively preserved. 

It is worth mentioning that for phenomenological boson models including high anhar- 
monicities exact results are lacking. In a previous publication we studied three fourth order 
boson Hamiltonians which are exactly solvable jl5|. For two of them we analyzed the possi- 
bility of classifying the states describing the intrinsic degrees of freedom in rotational bands. 

In the present paper we continue the project started in our previous paper and study semi- 
classically a sixth order boson Hamiltonian which is exactly solvable. Results for the exact 
trajectories, their periods and quantization are presented in an analytical form. Interesting 
results for the transmission probability through a barrier are obtained within the WKB 
approach. 
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The formalism and final results are presented according to the following plan. In Section 
II, a sextic boson Hamiltonian is chosen. This is treated within a time dependent variational 
principle approach in Section III. The quantization of one periodic coordinate and the semi- 
classical spectrum is given for several distinct energy regions in Section IV. The main results 
and final conclusions are collected in Section V. 



II. THE MODEL SEXTIC BOSON HAMILTONIAN 



In the present paper we shall treat semi-classically a sixth order quadrupole boson Hamil- 
tonian 

H^A^Y, b% + A2Y^ {blbl, + b,b,) {-r + A,P' + A,P\ (2.1) 

where (6^) with — 2 < < 2, denotes the creation (annihilation) quadrupole boson 
operator and 

^ = ^ E w + (^-m + ^-y^^) (-)"• (2-2) 

Writing the high order boson terms, involved in if, in a normal order it results that the 
model Hamiltonian does not commute with the boson number operator 

(2.3) 

In what follows we are interested in solving the time dependent variational principle (TDVP) 
equations 

S J^{^\(^H-ih^^\^f)dt' = 0. (2.4) 

When the variational state |^) spans the whole space of the boson states, solving the equa- 
tion (2.4) is equivalent to solving the time dependent equation associated to the model 
Hamiltonian H. The classical features encountered by H can be described by restricting the 
space of 1^) to the coherent states: 

I*) = exp [zobj - z*bo + Z2(b^ + hl^) - z*{h2 + b_2)] |0). (2.5) 

Here the boson vacuum state is denoted by |0). The function |^) depends on the complex 
parameters zo, Z2 and their complex conjugates Zq,Z2- These parameters play the role of 
classical phase space coordinates whose equations of motion are provided by the TDVP 
equations. 



It is convenient to have the classical equations of motion in a canonical form. In order to 
touch this goal we perform the following change of variables 

= 2^'^+'y'Re{z,), = 2('=+2)/^Im(z,), k = 0,2, i = (2.6) 

Up to an additive constant, ^^44+^^6, the classical energy function, defined as the average 
value of the model Hamiltonian on the coherent state |\E'), has the expression: 

H ^ {m\^) = ^{pI+pD + ^{ql + <&) + jiqf + qlf + f (g? + qI? , (2.7) 
where the following notations have been used: 
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A' = Ai- 2A2, A = A1 + A2 + UA4 + —Ae , 

D = AAi + 54^6, F = 6^6. (2.8) 

The nice feature of the classical energy function Ti. consists of the fact that it does not 
contain powers of momenta higher than two, although we started with a high order boson 
Hamiltonian. 

It is worth mentioning that similar time dependent approaches are also used for treating 
many body systems, when for example a Hartree Fock or an RPA approximation jl^ 0| is 
performed. In contradistinction to those formalisms where the classical coordinates, param- 
eterizing the trial wave function, are small expansion parameters, here no approximation is 
made in deriving the classical equations of motion from the variational principle of Eq.(2.4) 

The objective of the next section is to study the classical motion of the system associated 
to the classical energy function H. 



III. THE CLASSICAL DESCRIPTION 



The classical energy expressed in terms of the polar coordinates (r, 9) associated to the 
plane {qi,q2), looks like: 

n = ^{P + r'e') + V{r), (3.1) 

where V{r) is the potential energy which contains in addition to the potential used in the 
previous publication, a sextic term: 

V{r) = ^Ar^ + ^Dr^ + ^Fr^ . (3.2) 
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For r > 0, depending on the coefficients involved, the potential energy function exhibits 
either two or no extreme points. In the ffist situation the ordering of the maximum and 
minimum points is decided by the relative signs of the defining coefficients A, D and F. Here 
we study the case A>0, D<0, F>0 which defines a potential having ffist a maximum 
and then a minimum. The other ordering situation will be considered in a subsequent work. 
Throughout this paper the applications are made with the following values: 

A = 3 MeV, D = -0.4 MeV, F = 0.01 MeV, A' = 0.0025 MeV. (3.3) 
There are two constants of the motion, the energy: 

n = E, (3.4) 

and the third component of the angular momentum, 

h ■ 

^3 = -^{qm - q2Pi) = -^r^o ■ (3.5) 

For the sake of completeness we give also the expressions of the other two components of 
angular momentum 

^1 = ^{<ii+p\- ql-pl), ^2 = ^(qm+pm)- (3.6) 

These components generate a classical SUc{2) algebra with the multiplication operation: 

{Ci, Ck} = ^T'^ikj^j , (3.7) 

where {, } denotes the Poisson brackets and eikj the antisymmetric unit tensor. Ci might 
be obtained as average values on ( see Eq. (2.5)) of angular momentum components 
operators Lj acting on a boson space 15|. These operators generate a boson SUb{2) algebra. 
The correspondence [, ] —>■ realizes a homeomorphism of the two algebras mentioned 

above. The existence of the constant of motion is a consequence of the fact the classical 
equations of motion are obtained through the variational principle equation (2.4). On the 
other hand the equation £3 = const, holds since in the quantal picture the conservation 
law equation [H, L^] = takes place. The classical system has two constants of motion Ti 
and £3 and two degrees of freedom r and 6. Therefore the system is fully integrable and 
consequently analytical solutions for trajectories are expected. As a matter of fact this is 
one of the aims of the present investigation. 



Note that ^ plays the role of the mass , mo, of a classical non-relativistic particle moving 
in a central force field. Also we note that when the sixth power term in r is missing one 
obtains the classical energy function used in our previous publication [15]. Moreover the 
classical energy may be viewed as a counterpart of a microscopic Hamiltonian including a 
two body monopole-monopole interaction. In that case the coordinate r signifies the classical 
image of a collective microscopic coordinate jisl ]. 

Eliminating the angular variable in Eq. p.ip one obtains: 
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+ V,ff{r) = E, (3.8) 



2A 

where Veff (r) is the effective potential energy given by: 

Veff{r) = ^,+V{r). (3.9) 

Since £3 is the third component of the intrinsic angular momentum, this is the classical 
counterpart of the K quantum number of the liquid drop model. Having this feature in 
mind we quantize it by the restriction 

C3 = 2Lh, L = 0,1,2,... (3.10) 

Replacing £3 by its quantized expression, the effective potential becomes: 

V^ffmr) = —^ + V{r). (3.11) 

For L < 59, Veff{L]r) exhibits a two minima shape. For the critical value L = 59, the 
minimum of the left well gets unified with the maximum of Veff which results in arising an 
inflection point. For L = 25 and the set of parameters specified by Eq.(3.3), Vgff (L; r) was 
plotted as function of r in Fig. 1. 

From Eqs. (3.8), (3.5) and (3.10), one obtains the equations of motion for x = r"^ and 6, 
respectively: 



h 3 dx 

dt 



2 V A'F ^P,ix-E,Ly 



de = L\l J L - , (3.12) 



where 



7 



with 



P^ix; E, L) = -x'^ - ax^ - - 7 {E) x - 5{L) 



(3.13) 



?)D ?)A 
a = —<0, f3 = jr>0, (3.14) 

1{E) = -jr<0, S{L) = 12-^ > 0. 

For any L < 59 there are three energy domains defined as follows (see figure 1 for the 
case L — 25): 

a) UNDER, for E e [Veff min2(L), Veff minl(L)], 

b) BETWEEN, for E G [Ves ,„ini(L), Veg ^^.(L)], (3.15) 

c) OVER, for E e [Veif ^ax(L), +00). 

The labels for the three intervals will be hereafter abbreviated by U, B and O respectively 
whenever one wants to mention the fact that a given observable characterizes a certain energy 
region. The system motion is allowed only for the values of "x" where the polynomial P 
acquires positive values. Such allowed intervals are depicted separately for each of the energy 
regions U,B and O. 

A. Results for an "U" energy 

In the energy interval U, the motion is possible only in the right well (see Fig. 1). Indeed, 
inside U, the equation 

Pi{x;E,L) ^0, (3.16) 

has two real 

Xi = a{E,L), X2 = b{E,L), xi > X2, (3-17) 
and two complex conjugate solutions: 

X3^u{E,L)+iv{E,L)^xl, (3.18) 

It is conspicuous now that P4{x; E,L) >0 for x e [x2, xi]. 
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For what follows it is useful to introduce the equation 

v^X"^ + {ah - (a + h)u + {u^ + v'^))\ - ^(a - hf = 0. (3.19) 

whose solutions are conventionally denoted by \i 2{E, L) (A2 > Ai). We notice that the first 
equation ()3.12|1 provides the dependence of time on the radial coordinate r. The analytical 
result for this dependence is given in Appendix A. By an inversion operation one obtains r 
as a function of time. Inserting this in Eq.(3.5) and integrating the resulting equation one 
obtains 9 as function of time. The two functions of time r and 6 are periodic. However, their 
periods are different from each other. Moreover they are not commensurable which results 
in having open trajectories. Here we give details about the description of the motion of the 
variable r, the calculations for the 9 variable being performed by following similar steps. 

For r, the period of the motion, Tu{E,L), i.e. twice the time elapsed between two 
successive passages through the turnin g po ints situated at r^-m = ^/b{E, L) and r^ax = 
^ya{E, L), is given by the equation [l9, 20, 21 1 



ME, L) = .fi^/ ,f 1 . 1 ; 1; i)), (3.20) 

where 

Ac/(^, L) = {ab - (a + b)u + {u^ + v^))^ + (a - 6) V, (3.21) 

with a, b and p depending on the energy E and the quantum number L as shown in Eqs. 
(3.19) and (3.20). 2Fi{a,b; c; z) is the Gauss hypergeometric function and 

The classical action corresponding to the degree of freedom r, written in units of 27r/j, is 
given by the following integral: 



1 



E 



Iu{E, L) = -- Tu{w, L)dw. (3.23) 

^^"^J Veff min2{L) 

B. Results for a "B" energy 

Similarly, in the case B, the zeros of the polynomial P4{x;E,L), are positive numbers 
denoted by 

xi = a{E,L), X2 = b{E,L), (3.24) 
X3 = c{E,L), X4 = d{E,L), 



and ordered as follows 

xi > X2 > X3 > X4. (3.25) 

The solution for t as a function of r is given in Appendix A. The periods in both wells, 
left and right, are equal to each other. Indeed, one can prove that the two periods have a 
common expression: 



^ ' V A'F ^{a{E, L) - c{E, L)){b{E, L) - d{E, L)) ^2' 2' 

(3.26) 

where 

;.2,z. r,_ W,L)-h{E,L)){c{E,L)-d{E,L)) 
"^^""^ - {a{E, L) - c{E, imE, L) - d{E, L)) ' ^'-''^ 
Since the hypergeometric function has a simple pole in A; = 1 the period of the motion 
diverges when the energy approaches Veff max{L) from below. If we calculate the period for 
E — Veff maxiL), the Same divergence is obtained. 

For the left well, the integral action given in units of 27rh is equal to that formulated for 
the right well but restricted to the corresponding energy interval: 

1 /-^ 

7,eft(£;, L) = Ib{E, L) = -— Tb{w, L)dw. (3.28) 

The action for a trajectory of energy E lying in the right well consists of two terms: 

/nght(£^, L) = Iu{Veff minl(iv), L) + Ib{E, L). (3.29) 



It turns out that luiVeff mini(-f'), L) is an integer number, namely 

"Ve// minlW 
Veff min2(i) 

Thus , the final expression for the integral action for a trajectory from the right well is 



1 /-^e// minlW 

Iu{Veff rainliL), L) = -— T^{w, L)dw = L. (3.30) 



/right (£;, L) = L + 7ieft(£;, L). (3.31) 
We notice that k%{E,L) < 1 for < K// max(-^^), but k\{Veff rmx{L), L) = 1. 

C. Results for an "O" energy 

The case c) called " OVER" is similar to the case " UNDER" . The only change to be done 
in the formulae pertaining to case U is the mere replacement b{E, L) —>■ d{E, L) required 
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by the conventional designations of the zeros of the polynomial P4. Thus, the roots of 
P4(x; E, L), are 

xi = a{E, L), X2= u{E, L) + iv{E, L) = x^, X4 = d{E, L), (3.32) 

with xi > X4. (See also the notations for the case b) BETWEEN, when all roots are real 
and positive.) 

Correspondingly, we denote by \i^2{E,L) (A2 > Ai) the solutions of the equation 

v^X'^ + (ad - (a + d)u + + v^)X -^{a-df = 0. (3.33) 

For a chosen energy in the interval O, the first equation (j3.12p has the analytical solution 
given in Appendix A. 

The period of the motion of the coordinate r, to{E,L), i.e. twice the time elapsed 



between two successive passages through the turning points situated at r^^^ = ^Jd{E, L) 



and Tmax = V o(-E, L), is given by the equation 



where 



To(E, L) = .ft^/ J_^=l= 1 L l; kUE. L)), (3.34) 

and Ao{E,L) is the discriminant of Eq. ()3.33p . When the energy approaches K//max(-^) 
from above, A2 00, and kp 1. Thus, the period of the motion diverges also when the 
energy approaches Ve/f ma.x{L) from above. 
The integral action for the case O is: 

1 

Io{E, L) = /right(K// max(i^), L) + —- To{w, L)dw. (3.36) 



D. The virtual motion under the hump 

The region shown in Fig.l under the hump is forbidden for the classical motion of r. 
Indeed, there the energy is smaller than Vef] which results in having a negative kinetic 
energy as required by Eq. (j3.8|) . Note that for this situation, the integral (j3.12p provides 
imaginary values for time, i.e, r = it. It is interesting to notice that by changing the time 
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variable t by —ir the resulting equation can be again integrated since the corresponding 
effective potential is now — Vg// while the energy becomes -E. Trajectories in the well —Veff 
are conventionally called virtual. They are periodic curves in the new variable r. Moreover 
this period may be connected with the tunneling time through the hump, for the associated 
quantum system 22) 2^- 

The period of the motion in the well —Veff {L; r), or in the hump of Veff {L; r), is given 

by 



3 1 11 

T^irtiE,L) = TTh\l- — 2^i{-, -; l;kl(E,L)), 

(3.37) 

where 

, 2 n - iKE,L)-ciE,L)){aiE,L)-d{E,L)) 

" {a{E, L) - c{E, L)){b{E, L) - d{E, L)) ' ^^'^^^ 
The integral action for the inverse of the hump potential is obtained in an analogous way 

with the procedure described before. The final result is: 

\ r^eff max(i) 

/virt(^, L) = ^J ^ T^iAw, L)dw. (3.39) 

E. Numerical application 

In order to emphasize on the complex structure of the classical motion here we give some 
relevant results for trajectories and their periods. The coefficients used in our calculations 
are listed in Eq.(3.3). The shape of trajectories depend on the chosen value of energy. Here 
we consider a pair of trajectories corresponding to an energy lying close to the maximum 
value of the effective potential, Fig. 2. One trajectory is lying in the left well, while the 
other one in the right well. The two trajectories are quite different. Indeed, as shown in 
Fig. 2, in the left well r is running most of the distance ^/xl — ^Jx^ in a very short time, 
than spends a large interval of time to reach the value ^/x^ and to depart from it . Finally 
r is coming quickly back to the value ^/xi. The angle 9 undertakes a jump of 2.5 rad before 
ending a period of motion. As for the right well the r curve corresponding to E close to 
Veff max IS different from the one commented above. Indeed, Fig. 2 shows that r stays longer 
close to the value ^/x2 and a very short time in the neighborhood of ^Jxl. On the other 
hand, 9 changes the speed around the half of period which is at variance with the behavior 
of the trajectory from the left well. 
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From Fig. 2 we see that after an interval of two periods of r, 2T, the variable 9 covers 
about 11 rad in the left well and only 5 rad in the right well. This suggests that the trajectory 
r = r{9) is not closed. The two trajectories, lying in the left and right well respectively, 
correspond to E= 6.984326 MeV and are presented in Fig. 3. As shown there, the inner 
and outer trajectories are almost tangent to each other. When the energy value is far from 
Veff,max there is an unreachable region, due to the bump. The periods of trajectories lying 
close to the maximum value of Veff are much larger than those staying far from the top of 
the hump. As a matter of fact, the period of r has a discontinuity for '^=Veff maxi which is 
pictorially shown in Fig. 4. Such a discontinue behavior suggests that Veff max is a critical 
point of a phase transition. In one phase, for a given energy, the system may follow one of 
two r trajectories of equal periods, depending on the initial conditions, while in the second 
phase only one trajectory is possible. 

IV. THE QUANTIZATION OF THE PERIODIC TRAJECTORIES 

We note that the classical Hamiltonian given by Eq. ()3.1|) does not depend on the angle 
9 but only on its conjugate momentum. Consequently there is a constant of motion £3 
which is discretetized by the constraint ()3.10p which in fact is a way of quantizing it. The 
remaining variable r and its conjugate momentum pr are both involved in the classical 
Hamiltonian. Moreover, as we have seen in the previous Section, r performs a periodical 
motion. Consequently, its motion can be quantized by a constraint for the classical action 
similar to the Bohr-Sommerfeld condition. This condition will be applied separately for each 
region specified in Fig.l. Thus, in the energy region called "UNDER", the quantization 
equation 



gives the energy levels EunAciiL^ni) situated in the right well below Veff mini{L). It is re- 
markable that the above equation is satisfied by 



This equation says that Veff mini(-f') is a quantal level in both left and right wells. Moreover, 
below this energy, in the right well there are another L levels. The quantization condition 



Iv{E,L) = m 



(4.1) 



^ = ^e// minl(-^^) , = L. 



(4.2) 
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in the left and right wells of the region "BETWEEN", reads 



Iieh{E, L) = ^2, /right L) ^ni, (4.3) 

These equations determine the same energy levels in both wells 

Eieft(L, na) = Eright{L, ni), (4.4) 

where rii — n2 + L. The spectrum in the region " OVER" is obtained by solving the equation 
provided by the quantization condition 

Io{E,L)^ns. (4.5) 

The virtual quantum states lying in the potential -Veff are obtained from 

7,irt(^, L) = 714. (4.6) 

Having the integral action for the forbidden energy region, one can calculate the transmission 
coefficient through the potential barrier. Indeed, according to the WKB approximation, the 
transmission coefficient through the hump is given by the equation: 

D{E, L) « exp(-27r7virt(£;, L)). (4.7) 

A. Numerical results 

Solving Eqs. (4.1), (4.3) and (4.5) satisfied by E, one finds the semiclassical spectrum of 
the system. It is instructive to see how many quantum states with energy less than V^ffmax 
accommodates the left well. The number of states is presented in Fig. 5 as a function of 
angular momentum. Note that for some L, the number of states is the same as for the L+1 
case. The quantal states are tentatively organized in rotational bands in Fig. 6, in a full 
analogy with the standard definition of the lowest bands within the liquid drop model. As 
for the right well the number of states and the rotational bands are shown in Figs. 7 and 
8 respectively. Figure 7 shows that for L < 10 the number of states of energy less then 
^eff max IS 48, whilc for the next 12 angular momenta the number of states in the interval 
mentioned above is 49. Intervals for angular momenta for which the same number of states 
is registered is decreasing and contains only one angular momentum when the maximum for 
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the effective potential disappears. The "bands" in the right well differ the ones from the left 
well. Indeed, as shown in Fig. 8, the energy spacing is much smaller then in the left well. 
This is caused by the fact that the secondary minima of Vgff are very close to each other, 
which is consistent with the known fact that superdeformed nuclei exhibit larger moments 
of inertia. Another difference between the two sets of bands consists of that while in the left 
well the n = 2 band is lower in energy than the n = 1 band in the right well the ordering of 
the two bands is opposite. 

In contradistinction to what happens for real trajectories, the integral action for the 
virtual states, is a decreasing function of energy. Also, the number of states from —Vejf is 
a decreasing function of angular momentum (see Fig. 9). 

Energy levels for the virtual states do not coincide with energy levels corresponding to 
real states of left and right well. The differences in energy and the period of virtual states 
satisfy the uncertainty relations. 

Let us consider the transmission coefficients corresponding to the virtual state of highest 
energy, the first one under the hump. These are calculated within the WKB approximation 
by making use of Eq.(47), for a given value of the angular momentum. The results are 
given in Fig. 10. It is worth to mention the discontinuity of D(E,L) given by Eq. ()4.7|) 
for some values of L. Indeed, going from the values 10, 23, 29, 35, 40, 44, 47, 50, 53 to 
11, 24, 30, 36, 41, 45, 48, 51, 54, respectively, one achieves a very big jump in magnitude 
for the transmission coefficient. This result might have an important consequence for the 
mechanism of alpha or heavy cluster emission. Indeed, if the alpha's or heavy fragment 
would be in a state of angular momentum staying in the list from above, its penetration 
probability would be enlarged and a sizable emission rate could be registered. Having in 
mind the integration limits in Eq.(4.7), the transmission coefficient D(E,L), associated to 
the first virtual state depends on the difference between K//max and the state energy. The 
smaller the difference the larger the transmission coefficient. Increasing the difference quoted 
above, the transmission coefficient decreases more or less exponentially. Such a dependence 
may be seen in Fig. 11. 
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V. COMPARISON WITH THE PREVIOUS WORK. 



Although the present and previous work [15| use the same formahsm to treat two distinct 
Hamiltonians, the results are both technically and phy sically different from each other. 



Thus, the fourth order boson Hamiltonian used in Ref. [15| yields a classical effective potential 
which misses the superdeformed minimum. Consequently, while here one distinguishes three 
distinct energy regions, called O, B and U, with topologically different trajectories, in the 
quoted reference there is only one region of finite trajectories and that correspond to the left 
well from Fig. 2. Here all trajectories r(t) are periodical. The period depends on energy. 
Indeed, increasing the energy from the U to the B region in the right well, there is a smooth 
variation for period. The periods associated to trajectories of the same energy but lying in 
the left and right well respectively, are equal to each other. At the border between the B and 
the O regions the period is singular, which suggests a phase transition. In the case of four 
order boson Hamiltonian such a transition appear to take place between a phase of periodic 
orbits and one of open orbits. Although the periods for orbits of left and right B trajectories 
of similar energies are equal to each other, the time dependence of r and 9 coordinates are 
different, respectively. The difference is increasing when the energy approaches K//max- 
As a matter of fact this reflects the influence of the deformation on trajectories. In this 
context one may say that in ref. Q] only moderately deformed trajectories show up which 
results in having in the plot corresponding to Fig. 3 only the inner curves. Thus, the 
superdeformed trajectories are specific to the present formalism. The absence of periodic 
orbits for E > V^ffmax has as consequence the fact that in a plot similar to that from Fig. 
4, the descending branch does not exist. 

Since for the chosen parameters the energy levels are almost equidistant one may expand 
the effective potentials around the minima and keep only the quadratic terms. The frequen- 
cies characterizing the harmonic motion in the two resulting wells have the expressions: 



1/2 

, fc=l,2 (5.1) 



where r^ i and ro,2 denote the coordinates of the first and second minimum, respectively. 
Note that for k=l, the influence of the sixth order boson term is the Tq term in the above 
expression while the k = 2 case exists only if the sixth order term is included. For the two 
wells one may define rotational bands in a similar way as in the framework of the liquid 
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drop model. 

4°i = V:]f{L,ro,k). L = 0,2,4,...;A; = 1,2, 
4!1 = V;7}"(L,ro,,)+c^.,L, L = 2,3,4,...;fc = l,2, 



= KF(^'^o.'^) + 2^A.,L, L = 0,2,4,...;fc = l,2. (5.2) 

It is worth mentioning that the first equation ()5.2j) with k=l represents a generahzation 
of the variable moment inertia formula 0, |2^ . When Vg/ / is truncated to the quadratic 
term, the equation expressing the energy in terms of angular momentum derived in Ref. 2^ 
is obtained. The above cited equation considered for = 2 represents an extension of the 
variable moment of inertia formula to the superdeformed nuclei. The second and third 
equation from ()5.2|) describe wobbling motions around the states from the ground band. 



3 



since there, no right well in the 



The results shown in Figs. 7, 8 are missing in ref. 
effective potential shows up. Here we presented also the results characterizing the virtual 
states lying under the bump. Also the transmission coefficients through the bump for a 
particular energy has been studied as function of both angular momentum and the energy 
normalized to its maximum value. 

The peculiar feature of the model Hamiltonian studied so far consists of that it contains 
only a quadratic dependence on momenta which is consistent with the classical image where 
the kinetic energy involves only such type of terms. Moreover, six is the maximal order of 
boson terms which can be analytically treated. From this point of view the present treatment 
cannot be improved by adding higher order terms. 

Formalisms dealing with anharmonic boson terms use Hamiltonians having, however, 
powers of momenta higher than two. Two examples of this type have been considered 
in Ref. [isj] . Despite their complex structure, the chosen boson Hamiltonians were exactly 
solvable, producing finite formulae for energy as function of various quantum numbers. For 
the sake of saving the space here we consider only one of the two Hamiltonians mentioned 
above, corrected with a sixth order boson term: 

H2 = tY,b%+ J2 Cj[ib^b^)jibb)j)]^ + Fib^b%N{bb)o. (5.3) 

M J=0,2,4 

Here N denotes the quadrupole boson number operator (2.3). Using the same notations as 
in Ref. [iS^ . the average of H2 on |\1') (2.5) has the expression: 

^2 = ^(g? + ql +pI+ pI) + j{qf + ql + P? + pIY + j{qiP2 - q2Pif 
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{ql + ql+pl+pD- (5.4) 



As shown in R,ef . [ifij] . and £3 are constants of motion. Taking the two constants equal 
to h^L{L + 1) and HM respectively, the energy becomes a function of the quantum numbers 
L and M: 

Elm = 2A^L{L + 1) + ABL{L + 1) + ^M^ 

+ F^VL(L + 1) + 1) - M^] . (5.5) 

This energy is associated to the motion of the intrinsic degrees of freedom. Supposing that 
these degrees of freedom are only weakly coupled to the Euler angles, then the total energy 
may be written as: 

EjLM = Elm + SJ{J + 1) (5.6) 

where J denotes the angular momentum in the laboratory frame. On the other hand, 
the ground state energies can be obtained by averaging H on the states {\NvaJM)} with 
N = v = j. The same energies are obtained if in Eq. ()5.5|1 one substitutes: 

2^/L{L + l)^^, M = (5.7) 

Thus, energies of the ground band are given by the following equation: 

Ej = A- + B— + 5J(J + 1) + - — 
^ 2 4 ^ ^58 

= AiJ + BiJ"^ + CiJ\ (5.8) 

The three parameters formula is expected to describe a larger number of experimental ground 
band energies than the two parameter formula from Ref. 15| . To give an example we present, 
in Table 1, the results of a least square fit of the ground band energies for ^^^Th. From Table 
1 one notices that the results provided by Eq. (|5.8|) are in much better agreement with the 
experimental data than those of Ref. [15]. 

Alternatively, the quantal energy for the intrinsic motion can be obtained by quantizing 
the anharmonic plane oscillator and the third component of angular momentum. In this 
case the result is: 

En,M = A{n + 1) + B{n + 1)^ + ^M^ + ^ [{n + if - 4(n + 1)M^] . (5.9) 
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Ex. 


present 


Ref.[15] 




Ex. 


present 


Ref. [15] 


2+ 


0.049 


0.059 


0.106 


16+ 


1.859 


1.849 


1.823 


4+ 


0.162 


0.175 


0.246 


18+ 


2.263 


2.257 


2.208 


6+ 


0.333 


0.343 


0.421 


20+ 


2.692 


2.693 


2.629 


8+ 


0.557 


0.561 


0.632 


22+ 


3.144 


3.152 


3.084 


10+ 


0.827 


0.824 


0.877 


24+ 


3.620 


3.631 


3.574 


12+ 


1.137 


1.129 


1.158 


26+ 


4.116 


4.126 


4.099 


14+ 


1.483 


1.472 


1.473 


28+ 


4.632 


4.633 


4.659 










30+ 


5.162 


5.150 


5.254 



TABLE I: Excitation energies obtained with Ea. H5.8|) and with the two parameters formula derived 
in Ref. [15], respectively, and given in units of MeV, are compared with the corresponding experi- 



mental data 26] for ^^^Th. The coefficients are determined by a least square procedure. The results 
for the coefficients involved in Eq^^ are:Ai = 14.96872 keV A2 = 7.50384A;ey and ^3 = -0.07602 
keV while the two parameters from the corresponding equation of Ref. [15] aie-.Ai = 44.01753 keV, 
A2 = 4.37055keV. 



Neglecting the zero point energy terms and taking M = \/L{L + 1) one obtains: 

EnL = A,n - Bin^ + CiL{L + 1) + :^ (n^ - AnL{L + 1)) 

5 



(5.10) 



with evident notations for the coefficients Ai,Bi,Ci. From the above equation one sees 
that the sixth order term brings two types of corrections, one depending on and one 
depending on the product nL{L + l). Equation ()5.10|) . involving four free parameters, seems 
to be very useful to describe in a realistic fashion the excitation energies of states having the 
same angular momentum. Extended calculations for ground band excitation energies with 
Eq. ()5.8j) as well as for excitation energies of a set of states of similar angular momentum 
with Eq. (j5.1Up are in progress. The results will be published elsewhere. 



VI. SUMMARY AND CONCLUSIONS 



In the previous sections a sixth order boson Hamiltonian was semi-classically studied. 
The minimal action principle applied to the model Hamiltonian and a variational state ]\E') 
of a coherent type, provides a set of equations for the complex parameters zq, Z2, Zq, z^ which 
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play the role of classical phase space coordinates. In terms of polar coordinates {r,6), the 
classical equations of motion have a canonical form. There are two constants of motion and 
therefore the system is fully integrable. The reduced classical phase space is one dimensional, 
the independent degree of freedom being the radius r. For a suitable set of the Hamiltonian 
coefficients, the effective potential involved in the classical energy function has a two minima 
shape. The differential equation of r can easily be integrated. The trajectories r{t) are 
periodic. Fixing the constant of motion C3 the equation of motion can be also analytically 
solved. The corresponding solution 6{t) is also a periodical function of time. The periods of 
r{t) and 6{t) are different from each other. Moreover they are incommensurable and therefore 
the trajectory r{6) is not closed. The emphasis is put mainly on the r(t) trajectories. There 
are three distinct energy regions, labeled by U, B and O, where the r{t) trajectories have 
distinct behavior. Indeed, for a given energy E and a given value of L (half of £3) in the 
region U there is only one periodic trajectory, in B the are two trajectories of equal periods 
while in O again there exists only one periodic trajectory. The period magnitude depends 
on energy. It is an increasing function for E < Vej/max and is singular at the interval 
end, while for the complementary region the period decreases when the energy is increased. 
This feature suggests that E = Veffmax is a critical point where the system undertakes a 
phase transition. Although the trajectories r{t) have a complex structure which is far from 
that of a periodic harmonic function, they can be quantized by a constraint for the integral 
action similar to the Bohr-Sommerfeld quantization condition. In this way one obtains a 
semi-classical spectrum for each of the regions U, B and O. Choosing the constraint 
for £3 in the region U there are L + 1 states while in B an equal number of states in the left 
and right wells, respectively. The states in B are tentatively organized in rotational bands 
keeping close to the conventional definition adopted by the liquid drop model. The states 
of the right well are very much compressed compared with those of the left well. Moreover 
the n = 2 band from the left well is lower in energy than the n = 1 band while in the right 
band the bands ordering is opposite. 

The tunneling phenomenon from the left to the right well through the separating barrier 
was also studied. A special attention is paid to the first virtual state under the hump. It 
is remarkable the fact that there are several values for the angular momentum L where the 
transmission coefficients, evaluated within the WKB approach, exhibit a spectacular increase 
in magnitude. It is interesting to see whether this result has some relevance for alpha or 
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heavy cluster emission. Indeed, if the preformed fragment has an angular momentum equal 
to one of the values where the transmission coefficients are very large one expects that the 
emission rate is adequately increased. Represented as a function of the energy distance 
between the chosen state and V^ffmax , the transmission coefficient is an exponentially 
decreasing function. 

A detailed co,„pa..o„ of the ptese. and Ref.Q's results ,s .ade Kav,„g the a„„ of 

pointing out the effect caused by the sixth order term. In order to have a complete view 
of this comparison additional equations for the sixth order Hamiltonian have been derived: 
a)one giving the harmonic energies in the two wells of the effective potential b) a three 
parameters formula for the ground band excitation energies and c) an equation for energy 
describing the intrinsic motion, depending on a plane oscillator quantum number and an 
angular momentum associated to the boson degrees of freedom. 



VII. APPENDIX A 

Here we give the analytical solution of the first equation ()3.12j) for an energy value be- 
longing to one of the three domains defined in Section III. Aiming at a compact presentation 
of the final results it is useful to introduce the following notations: 



Si = hJ^ ^ -. = Sr. (A.l) 

Here A', F are two of the Hamiltonian coefficients. A^- the second rank equation ()3.19l3.33p 
discriminants and a,b,c,d are the roots of the fourth order polynomial P^ix; E, L) introduced 
in previous sections. Two particular values of the variable x are to be specified: 

Ai + i A2 + i 

Ai and A2 are the solutions of Eqs. ()3.19p and ()3.33|) if m is equal to U and O, respectively. 
Also the following intervals for r are needed 



= [\/rf,v^, i?g^ = [Av^, (A.3) 
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Integrating the first equation ()3.12j) one obtains the time t as a function of r. In the intervals 
with 171 = 11,0 and R^^ with i=l,2 the final result can be written in a compact form: 



tm{E,L;r) = S^F{^imiE, L;r), kn{E, L)), m = U,l,r,0, 

n = m{6m,U + Sm,o) + B{6m,l + Sm,r)- 

If r belongs to -Rm"* with m = U,0, the solution is: 

UE, L- r) = 2Fi{^, ^; 1; k^E, L)) + F(02m(i?, L- r), k^{E, L)) ) . 

F(<l>, k) denotes the incomplete elliptic integral of the first kind 



Vl- /t^sin^e' 



(A.4) 



(A.5) 



(A.6) 



The functions $im and $2m, m = U,0 have the expressions 



L; r) = arccos 



')2m{E,L]r) = arcsin 



Vm - / Ai(A2 + 1) 



(A2 + l)(Ai-A2) 1 



A2(Ai + l) 

For the B energy region the functions $ are 



v 



Ai+1 / 



m = U,0. (A.7) 



hi{E,L;r) 
'^,{E,L;r) 



arcsin \ 



arcsin \ 



(a — c)(r2 — d) 
(c — d){a — r^) 

I — c)(r2 — b) 
(a — 6)(r2 — c) 



(A.8) 
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FIG. 1: The effective potential given by H3.9|) and ()3.2() with the coefficients A',A,D,F specified 
in Eq. is plotted as function of the adimensional variable r for L=25 where L is defined by 

Eq. (L3). 
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400 800 1200 1600 2000 2400 2800 400 800 1200 1600 2000 2400 2800 



ltf't[s] 10^t[s] 

FIG. 2: The trajectories r(t)(adimensional) and ^(t)(in units of rad) for the left and right wells, 
E=6. 984326 MeV and L=25, where L has the meaning specified by Eq. (|3.1()|) . 
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FIG. 3: The trajectories r(^)(adimensional) are represented in the Cartesian coordinates x = 
rcos{0) and y = rsin(^) for L=25 and E=6. 984326 MeV. The curves in both the left (the inner 
curve) and right (the outer curve) wells are given. 
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FIG. 4: The trajectory periods in the left well of the region B as well as in the region O are plotted 
as function of energy. 
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FIG. 5: The number of states from the left weh are given as function of L p.lU|l . 
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FIG. 6: The energy levels of the n = 0, n = 1 and n = 2 bands from the left well. 
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FIG. 7: The number of states in the right well is plotted as function of L defined by Eq. p.lU|l . 
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FIG. 8: The energy levels in the n = 0, n = 1 and n = 2 bands obtained in the right well. 
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FIG. 9: Number of virtual states which appear under the effective potential barrier. 
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FIG. 10: The transmission coefficient calculated with Eq. (4.7) is plotted as function of L given 

byEq. (irrm . 
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FIG. 11: The transmission coefficient (adimensional) defined by Eq. (|4.7|) is plotted vs. 
AEfirst-max the excitation energy (given in units of keV) of the first virtual state. 
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